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Abstract 

The heavy-quark spectral function in a hot plasma is reconstructed from the corresponding eu- 
clidean propagator. The latter is evaluated through a path-integral simulation. A weak-coupling 
calculation is also performed, allowing to interpret the qualitative behavior of the spectral func- 
tion in terms of quite general physical processes. 



1. Introduction 

Heavy Quarks (HQs), being produced in the very early stages of nucleus-nucleus collisions, 
have long been used as probes of the resulting hot (possibly deconfined) medium. Medium 
modifications of heavy-quark spectra (elliptic flow, suppression of high-p^ particles) can shed 
light on the degree of their interaction and thermalization with the rest of the plasma. As a 
result, it is interesting to perform a first-principle calculation of the heavy-quark spectral function 
in a hot plasma, showing how for the latter an interesting qualitative structure arises from the 
interaction with the medium. We first perform a Hard Thermal Loop (HTL) weak-coupling 
calculation, allowing to identify quite general physical processes leading to a modification of the 
vacuum result. We then formulate the problem in terms of a path-integral which is evaluated 
exactly, with a Monte Carlo algorithm. From such an euclidean correlator, we can shed light on 
the structure of the corresponding non-perturbative in-medium spectral density. 

2. HQ spectral function: weak-coupling calculation 

The in-medium analytic propagator of a non-relativistic heavy-quark reads 

G(z,p)= ~* , (1) 

Z-E p -I.(z,p) 
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Figure 1: The different processes contributing to the imaginary part of the HQ self-energy. 



where E p = M + p 2 /2M, S is the HQ self-energy and setting z = to + ir] gives the retarded 
propagator G R (to) = G(to+i rj). The imaginary part of the latter provides the HQ spectral function 

T(to) 



o-(to) = 2lmG R (to) = 



[cj-E„- ReS(w)] 2 + T 2 (w)/4' 



(2) 



where r(to) = -2lml. R (to). The HQ spectral function is then non-vanishing only for energies 
for which the self-energy develops an imaginary-part. For a non-relativistic heavy quark at zero 
momentum (playing the role of an external probe) placed in the QGP the self-energy can be 
evaluated in the (HTL-resummed) one-loop approximation 



1 + N(k°) 



N(k°) 



(3) 



p° - (E k + k°) p o - ( Ek - jfcO) 
which is expressed in terms of the HTL spectral function of a longitudinal gluon [1] 

Pl (oj >0,q) = 2tt [z L {q)6{co-co L {q)) + 0(q 2 - a?) p L {io,q)] . (4) 

The latter displays a peak for time-like momenta, arising from the propagation of a plasma 
wave, and a continuum part for space-like momenta, describing the Landau damping (i.e. soft 
collisions with the plasma particle). The imaginary part of the HQ self-energy gets then a pole 
and a continuum contribution from the different processes shown in Fig. 1. The HQ can go 
on-shell through the emission/absorption (with obvious meaning of the statistical factors) of a 
plasmon 

/dk 
— (2n)Z L (k)x 

x { (1 + N(co L (k))) 6 [co- (E k + co L (k))] + N{u L {k)) 5[co- (E k - co L {k))} } (5) 
or of a space-like (resummed) gluon exchanged with the medium particles 

V*>»\oS) = g 2 C F f f k dk a p L (k°,k)x 

x (2k) {[l + N(k°)] 6 [to- (E k + k )] + N(k°) d[to- (E k - k )]) . (6) 

The results for the corresponding in-medium spectral function are displayed in Fig. 2. Its main 
features are: a broadening and a negative shift of the main peak and the appearance of secondary 
peaks at energies corresponding to a huge density of states for the absorption/emission of a 
plasmon. The bumps arise from divergences in rP° le (w) (Van-Hove singularities). 
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Figure 2: The imaginary part of the HQ self-energy in the QGP and the resulting spectral function for different tempera- 
tures, setting M= 1.5 GeV and a, = 0.3. 



3. HQ spectral function: path-integral approach 

In order to formulate the problem in terms of a path-integral we employ the following strat- 
egy, borrowed from [2, 3], and extended to the case of finite-mass quarks. We consider the 
propagation of the HQ in the background of a gauge field, summing over all its possible trajec- 
tories and integrating over the field configurations with an action accounting for the presence 
of an hot medium. The field configurations are thus weighted according to the HTL effective 
action, which allows to properly dress the propagation of soft (long- wavelength) modes with the 
interaction with the plasma particles. Neglecting non-abelian effects one takes the above action 
to be gaussian and this allows to perform the functional integral over the field configurations ex- 
actly, thus reducing by orders of magnitude the required numerical efforts with respect to a 1QCD 
calculation. In practice we are addressing the simpler case of a hot-QED plasma, for which the 
above procedure is exact, fixing at the end the parameters entering in the calculation (tempera- 
ture, coupling and HQ mass) to values of phenomenological interest for the study of the QGP. 
We make the choice of losing some peculiar feature of QCD, to investigate exactly (without any 
T/M expansion) and with a huge statistics (not conceivable in a 1QCD simulation) important and 
very general medium effects. 

We thus perform numerical simulations for the following euclidean correlator, describing the 
propagation of a HQ from (0, 0) to (t, r), 



G(r,r)= f ' [£>z]exp - f dr' (m +\-Mz 2 \ 
Jz(0)=o' L Jo V 2 ] 

[dr' [ 
Jo Jo 



x exp 



~ ,' ... , dr"A T L (r' - t'\z(t')-z(t")) 



(7) 



which results expressed in terms of the HTL longitudinal gauge-boson propagator (after subtrac- 
tion of the instantaneous vacuum Coulomb interaction) 

Al(t, q) = f ^e-^pdqo, q)[6(r) + N(q )], (8) 

with pi given in Eq. (4). By integrating over all possible final positions we project it to zero 
momentum and reconstruct the spectral function of a HQ at rest in the hot plasma, after inverting 

3 



CoilStLIIll 

Lorentzian 
Gaussian 




6 _ 
I 60- 



T=200 MeV 
7=250 MeV 
T=300MeV 
T=350 MeV 
T=400 MeV 



0.2 0.4 0.6 0.8 



Figure 3: The HQ spectral density from the MEM procedure, for various choices of the default model at T = 200 MeV 
(left) and for different temperatures with a gaussian default (right). We set M— 1.5 GeV and a s =0.3. 



the relation 



G(r,p = 0) 



/da> 
2^ 



e'^oico, p = 0). 



(9) 



For this we employ a Maximum Entropy Method [4] algorithm. The latter allows to find the 
most probable spectral density compatible with the data and some prior knowledge embodied in 
a default model. In Fig. 3 we display the results obtained with different choices of the default 
model and for various temperatures. With a flat default model one gets simply a very broad bump, 
with a peak slightly shifted to lower energies. On the contrary, inserting richer information into 
the prior knowledge - e.g. the width and mass shift arising from the static (M = oo) case and 
the request of fulfilling a few sum rules - the spectral density displays a qualitative behavior 
(with small bumps away from the main peak) reminiscent of the one found in the perturbative 
calculation. Remarkably, these findings arise from the inversion algorithm itself, due to the 
high quality of the data, without putting any information on the presence of secondary peaks in 
the default model. In particular the shift of the low-energy bump with the temperature appears 
consistent with its physical interpretation in terms of plasmon-absorption. 



4. Conclusions 

We displayed results for the HQ spectral function in a hot plasma, like the QGP. A HTL weak- 
coupling calculation allowed to identify important qualitative features of the latter (broadening 
and shift of the main peak and appearance of secondary bumps), relating them to well-defined 
physical processes. The problem was then addressed non-perturbatively, through a path-integral 
simulation. Also in this case, a non-trivial structure was found in the spectrum. Further details 
will be given elsewhere [5]. It will be of interest to compare our findings with 1QCD results [6]. 
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